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MULTI-BUMP GROUND STATES OF THE FRACTIONAL 
GIERER-MEINHARDT SYSTEM IN R 

JUNCHENG WEI AND WEN YANG 


Abstract. In this paper we study ground-states of the fractional Gierer- 
Meinhardt system on the line, namely the solutions of the problem 

( (-A ) s u + u - £ = 0, in R, 

(— A) s v + e 2s v — u 2 = 0, in R, 

u,v > 0, u, v —> 0 as \x\ —> +oo. 

We prove that given any positive integer k, there exists a solution to this prob¬ 
lem for s G [2,1) exhibiting exactly k bumps in its u— component, separated 

from each other at a distance 0(e 4 3 ) for s E (^,1) and O(| loge:| 2 ) for 
s = 2 respectively, whenever e is sufficiently small. These bumps resemble the 
shape of the unique solution of 

(—A) S U + U — U 2 = 0, 0 < U(y) -> 0 as \y\ —>• oo. 


1. Introduction 

In this paper we consider the following fractional Gierer-Meinhardt system in R 

f {-A) s u + u- =0 inR, ( , 

\ (—A) s v + e 2s v — u 2 = 0 in R, 

where (—A) s ,0 < s < 1, denotes the fractional Laplace operator. (For the defini¬ 
tion, See Section 2 below.) 

When s = 1, this is the classical Gierer-Meinhardt system proposed by Gierer- 
Meinhardt in m in 1972. More precisely they considered the following reaction- 
diffusion system as a model of biological pattern formation 

! a t = dAa — a + ^, inf2x(0,£), 

ht=DAh — h + a 2 , infix(0,f), (1.2) 

d„a = d„h = 0, on dfl x (0, t), 

where d, D > 0 are diffusion rates, f 1 C R" is a bounded domain and d v denotes 
the derivative in the outer normal direction. The Gierer-Meinhardt system was 
used in [12] to model head formation of Hydra , an animal of a few millimeters in 
length, made up of approximately 100,000 cells of about fifteen different types. It 
consists of a ’’head” region located at one end along its length. Typical experi¬ 
ments with hydra involve removing part of the ’’head” region transplanted area is 
sufficiently far from the (old) head. These observations led to the assumption of 
the existence of two chemical substances a slowly diffusing activator and a rapidly 
diffusing inhibitor, whose concentrations at the point x G and time t > 0 are 
represented, respectively, by the quantities a{x , t) and h(x , t). Their diffusion rates, 
given by the positive constants d and D are then assumed to be that d <C D. The 
Gierer-Meinhardt system falls within the framework of a theory proposed by Turing 
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m in 1952 as a mathematical model for the development of complex organisms 
from a single cell. He speculated that localized peaks in concentration of chemical 
substances, known as inducers or morphogens, could be responsible for a group of 
cells developing differently from the surrounding cells. Turing discovered through 
linear analysis that a large difference in relative size of diffusivities for activating 
and inhibiting substances carries instability of the homogeneous, constant steady 
state, thus leading to the presence of nontrivial, possibly stable stationary config¬ 
urations. Activator-inhibitor systems have been used widely in the mathematical 
theory of biological pattern formation mm- Substantial research concerning this 
system has been generated in recent years. We refer the reader to the survey papers 
[24] . [35] and the book [36 ; for the overview of the subject. 

In the last twenty years there have been intensive research on the existence and 
stability of steady states of the Gierer-Meinhardt system (H in a bounded domain. 
It is known that there are multiple spikes solutions which may be stable. We refer 
to papers [ 6 ], [IS] , [30], ]28] , [331 ■ El] and the book [36] and the references therein. 
In the case of the domain being the whole space, after suitable rescaling the steady 
state problem of in becomes 

( Au - u + ^ = 0, u > 0, in R n 

< Av — e 2 v + u 2 = 0, v > 0, in R n (1.3) 

[ v —y 0 as |rc| —» +oo. 

A solution to m is called ground state. In the real line case the existence of single 
and multiple pulse solutions is proved independently by Doelman-Gardner-Kaper 
in [7] (via geometric dynamical system method) and by Chen-del Pino-Kowalczyk 
in [3 (via PDE reduction method). Similar results have been obtained for the case 
n = 2 by del Pino-Kowalczyk-Wei in [S], In higher dimensional case another type 
of solutions exist: solutions which are radially symmetric but have rings concentra¬ 
tions. We refer to Ni-Wei [25], Kolokolnikov-Wei m and Kolokolnikov-Wei-Yang 
[20] , In R 3 there exists also axially symmetric solution with smoke ring concen¬ 
trations. See Kolokolnikov-Ren [23j. The presence of steady configurations in the 
whole space appears driven by smallness of the relative size e 2 = of the diffusion 
rates of the activating and inhibiting substances. The present paper deals with the 
case of equation (11.31) in the nonlocal diffusion-fractional laplacian case (11.11) for 
s € [ 5 ,1)- Next we briefly discuss the fractional laplacian and nonlocal diffusion. 

In probability, we consider the random walk for Levy processes: 

u " +1 = E p ^ 

k 

where Pjk denotes the transition function which has a tail (i.e., power decay with the 
distance |j — /c|). By taking the limit we get an operator (— A) s as the infinitesimal 
generator of a Levy process: if X t is the isotropic 2s—stable levy process we have 

(—A ) s u(x) = lim — E[rt(a:) — u(x + X/41. 

0 + h 

When s = 1 it corresponds to the Brownian motion. 

Fractional diffusion equations have been used to model anomalously slow or fast 
scattering of particles in a variety of natural applications. A consideration of the 
problem of anomalous sub-diffusion with reactions in terms of continuous-time ran¬ 
dom walks (CTRWs) with sources and sinks leads to a fractional activator-inhibitor 
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model with a fractional order temporal derivative operating on the spatial Lapla- 
cian. A similar type of system has also been proposed for diffusion with reactions 
on a fractal. The problem of anomalous super-diffusion with reactions has also been 
considered and in this case a fractional reaction-diffusion model has been proposed 
with the spatial Laplacian replaced by a spatial fractional differential operator. If 
the reaction time is not short compared with the diffusion time in sub-diffusive 
systems with reactions (for example, if many encounters between reactants are re¬ 
quired before reactions proceed), then an alternate model has been proposed where 
the fractional order temporal derivative operates on both the spatial Laplacian and 
the reaction term. An important distinction between the two anomalous reaction- 
diffusion models becomes apparent when the concentration of species is spatially 
homogeneous, this latter model does not reduce to the classical macroscopic rate 
equations except when the diffusion is also non-anomalous. For the more back¬ 
ground on fractional reaction-diffusion system, we refer the readers to mmm 
and references therein. We shall mention that weakly nonlinear analysis has been 
done to constant equilibriums for Turing’s system by Henry-Langlands-Wearne |14] 
and Golovin-Matkowsky-Volpert [13]. In [23] the author studied a slightly different 
Gierer-Meinhardt system with fractional diffusion 

( a t = -d(-A) s a - a + in (-1,1) x (0, t), 

< h t = DAh — h + a 2 , in (—1,1) x (0, t), (1.4) 

[ a x (±l,t) = h x (±l,t) = 0 

where (—A) s denotes the spectral fractional Laplacian. 

In this paper we shall consider the existence of nonlinear patterns for the classical 
Gierer-Meinhardt system in the real line m with fractioal diffusion. A similar 
notable feature of problem m (as the classical case) is, it will be shown in this 
work, the presence of a large number of solutions (modulo translations) as the 
parameter e gets smaller. More precisely, given any positive integer k we find a 
number e such that if 0 < e < £k then there exists a solution exhibiting exactly 
k bumps in the activator. Besides, after appropriate re-scaling of u, these bumps 

1 —2s 

are approaching a universal profile and are separately from each other e ~“ and 
(—loge )2 for the case s € (^,1) and s = \ respectively. We remark that this 
phenomenon is intrinsic to the full system, for only one ground state of the following 
equation modulo translations exists ([9]) 

{-A) S U + U-U 2 =0 inR, 0 < U(y) -)■ 0 as \y\ ->■ oo. (1.5) 

Before stating our main results on the existence of such solutions, we need the fol¬ 
lowing preparations. In the sequel by U{x) we denote the unique radially symmetric 
solution of (11.51) . (For the existence and uniqueness, we refer to Frank-Lenzmann 
[9].) In order to describe the position of bumps, we set 

m 

3 s(5i, 52, • ■ ■ ,5m) = ^2(a s U( 2 9i )+ /? s £ 25-1 |25 l | 2s - 1 ) 
i =1 

+ 5Z ( a sU(q l - 5j) + /3 s £ 2s_1 |5i - 5 j| 2s_1 ) 

+ (a s U{qi+qj) +/3 s e 2s ~ 1 \q i +q j \ 2s ~ 1 ^j 


(1.6) 
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for s € ( 2 ; 1 ) and 

m ^ 

Si(gi,g 2 ,--- ,9m) = (^1^(291) + / 3 i -—y-log| 2 gj|) 


£=1 




+ 51 ( a i ^^ - 9j) + /3i !°g l9i - Qj l) 

1 


+ Y1 + &) +/ ? i^T 1 °g|9i + 9j|) (1-7) 




for s = |. 

In Section 8 , we shall prove the function S s , s € [|,1) admits global minimal 
point in the interior of the following set 

Qs, V = {( 91 ,92,"' ,9m) I > 9* > T )£ L ^ S ', |9i - 9j| > 776^1 (1.8) 

for s € (|, 1 ) and 


Q±,r, = {( 91 ,92," • ,9m) | i(logi)= > Qi > 7?(log i)5, | ?i - ^ | > ^(log^j 


(1.9) 


for s = 1 respectively. Here 77 denotes a small positive number and with the 
constants a s ,/3 s being given later. 

We denote one of the global minimal points (if there are more) of 5 s in Qs, v by 

Qs (9s,1,9s,2, , 9s,m)' 

Remark: Since the function E s (q^ , 92 , • • • ,9m) is analytic in Q s v , therefore, all 
the global minimal points are discrete. 


Let us set 

r= j ( 2 ssin(jc) gl ~ 2s Jr 1 U 2 {y)dy^ , s G (|, 1), 
[ (^ lo g|/Ri^ 2 (»)dtf) > s = 


Our main result is the following. 


Theorem 1.1. Let k = 2m, m > 1 be a fixed positive integer. There exists £k > 0 
such that for each 0 < e < £*,, problem EH admits a solution (u, v) with the 
following property: 


lim 
£—>-0 


r c l u< 


■{x) -^2(U(x- q s ,i) +U(x + q Sti )) 


= 0, 


uniformly in x € M, while for the second component v, we have 


lim ( r e 1 v e (q s . i + x) - 1 

E-X) V 


-1 


for any i uniformly in compact sets in x, i = 1,2 
when k = 2m + 1 , m > 0 . 


v E (q s ,i - a;) - 1 ) =0 

,m. A similar result holds 
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We will only give the details of constructing solutions exhibiting an even number 
of bumps in the line. The case of odd number of bumps can be treated in a similar 
manner, and we will discuss the necessary changes in section 8 . 


The method employed in the proof of Theorem 11.11 consists of a Lyapunov- 
Schmidt type reduction. Fixing m points which from Q s ,ri with ry being determined 
later, an auxiliary problem is solved uniquely, and solutions satisfying the required 
conditions will be those precisely satisfying a nonlinear system of equations of the 
form 

c s ,i(gi,<? 2 , • • • ,9m) =0, i = 1,2, • • • ,m, 
where for such a class of points the function c Sj i satisfy 


g r i 

Cs,i (?1 , 92 , ' ' ' , 9m) = " 2 ^>(1^*1) ^ ( (-^s (1 9* 9.7 I ) "b F s (I <7,; + Qj |) ) 

1 

where 


u s,n 

(i.ii) 


and 


a s U (r) + f3 s £ 2s 1 r 2s x , s 6 ( 5 , 1 ), 
aiU(r) + P^, s = l 


( o(l)e s+ 2 2s , se(±,l), 

1 o(l)(-loge)"i, s = 5 - 


We can easily find that solutions of the problem c s> j = 0 are closely related to 
critical points of the functional defined in (II.61) and for s € (^, 1 ) and s = ^ 
respectively. In the classical case, since the profile U of each bump is exponential 
decay, in the leading term of Ci, we only need to consider the neighbor points of 97 . 
However, in the fractional case, the profile U has algebraic decay. As a result, all the 
points interact with each other strongly and we have to consider all the interactions 
in dealing with c s y. This is the main difficulty and new feature in considering cu¬ 


lt is clear that the results of Theorem 11.1 1 can be extended without any difficulty 
to Gierer-Meinhardt system in which activator and inhibitor have different diffusion 
characters 

f {~A) Sl u + u- ^ =0 inR, m 

\ (— A) S2 v + £ 2s2 v — u 2 = 0 inR, 

where si € (0,1] and s 2 G [|, 1]. It is also possible to generalize to Gierer-Meinhardt 
system with more general nonlinearity. We omit the details. 


The rest of the paper will be devoted to the proof of Theorem ll.il In section 2 we 
study the fractional Laplacian operator and the behavior of the Green function for 
(—A) s + 1. In section 3 we set up the scheme of proof, in particular we explain why 
the constant r e is the right scaling factor to get the desired multi-bump expansion. 
The program outlined there is carried over the following sections. 
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2. Preliminaries 


In this section, we provide some elementary properties of the operators (—A ) S +I. 
Let 0 < s < 1. Various definitions of the fractional Laplacian (—A) s </> of a function 
<j> defined in R" are available, depending on its regularity and growth properties. 

For <f> £ H 2s (M. n ), the standard definition is given via Fourier transform ~ . 
(—A) s (^> £ L 2 (R") is defined by the formula 




( 2 . 1 ) 


When (f> is assumed in addition sufficiently regular, we obtain the direct represen¬ 
tation 

(-a = 

J K" 


-d y 


( 2 . 2 ) 


|a: - y\ n + 2s 

for a suitable constant d SyU and the integral is understood in a principal value sense. 
This integral makes sense directly when s < ^ and (j> £ (7°’“(R") with a > 2s, or if 
<f> £ C' 1,Q (R n ), 1 + a > 2s. In the latter case, we can desingularize the integral and 
represent it in the form 

<t>(x) - 4>{y) ~ V0(x)(x - y) 


(—A ) s (j)(x) = d Stn f 

J R r 


\n-\-2s 


-d y. 


/i” \x - y\ r 

Another useful (local) representation, found by Caffarelli and Silverstre [2], is via 
the following boundary value problem in the half space R" +1 = {(x, y) \ x £ R n , y > 
0 }: 

V • (y 1-2s V(/>) = 0 in R^ +1 , 

4>(x, 0) = cj>(x) on R”. 

Here <j) is the s—harmonic extension of q i, explicitly given as a convolution integral 
with the s—Poisson kernel p s (x, y ), 


4>{x,y) = Ps( x - z,y)<j)(z)dz, 


where 


p a (x,y) = C n , s 


,,4s—1 


(N 2 + M 2 


and C n . s achieves J Rn p(x,y) = 1. Then under suitable regularity, (—A ) s (j) is the 
Dirichlet-to-Neumann map for this problem, namely 

(-A ) S( />{x)= lim y 1 ~ 2s d y (f>(x, y). (2.3) 

3/-> o+ 

Characterizations (HOT)-(12751) are all equivalent for instance in Schwartz’s space of 
rapidly decreasing smooth functions. 

Now let us consider for a number m > 0 and g £ L 2 { R") the equation 

(—A) s </> + mcj) = g in R ra . 

Then in terms of Fourier transform, this problem, for cf> £ L 2 , reads 

(l£l 2s +m)4> = g 

and has a unique solution <fi £ 7L 2s (R n ) given by the convolution 


4>{ x ) = T m[g ] ~ [ G(x - z)g(z)dz, 

JR n 


(2.4) 
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where 

^ = |£| 2s +m' 

Using the characterization (12.31) written in weak form, <f> can be characterized by 
(f>(x) = 0) in trace sense, where <j> G H is the unique solution of 

/ / V0V(^y 1_2s + m / tjxp = I gip, for all tp G H, (2.5) 

J JR" +1 JI” J R» 

where H is the Hilbert space of functions ip G H 1 ) )C (1R" +1 ) such that 

IMI h~ [ [ IV^IV 2s +m/ M 2 <+oo, 

J JR" +1 J R" 

or equivalently the closure of the set of all functions in G£°(M? + ) under this norm. 

A useful fact for our purpose is the equivalence of the representations flUD and 
(12.51) for g G L 2 (R n ). 

Lemma 2.1. Let g G L 2 (W rl ). Then the unique solution <fi G H of problem I A 51) is 
given by the s—harmonic extension of the function <f> = T m [g\ = G * g. 

For a proof, one can see Lemma 2.1 in [4j- Let us recall the main properties 
of the fundamental solution G(x) in the representation (12.41) . which are stated for 
instance in m and uni. 

We have that G is radially symmetric and positive, G G C°°(M n \ {0}) satisfying 
l<3(x)| + |x||VG(x)| < f_ 2a for all |x| < 1, 
lim G(x)|x| n+2s = 7 > 0, 

| x | —yoo 

|x||VG(x)| < ^ 2s for all |x| > 1. 

In next lemma, we get a more specific behavior of the Green function G(x) 
around x = 0 for n = 1 . 


Lemma 2.2. Let G(x) be the Green function of the following equation 

(—A) s G(x) + G(x) = 8(x). 

Then, we have 

qq + ai|x| 2s_1 + 0(|x| mm i 2 ’ 4s_1 }) as |x| —> 0, if s G (|, 1), 


G(x) = 


if s = 


lo g |*| + a 2 + 0(|x|) as |x| 0 

where ao, ai, 02 will be given in the proof. 

Remark: cii is negative. 

Proof. By using fourier transform, we can write the equation (El as 

(|^ + 1)G(0 = 1. 

Therefore, we have 


G(x) = 


1 


0 ix(, 


2tT 7-00 1 + l£| 2s 77 


cos(x£) 

T+TsF 


dG 


( 2 . 6 ) 

(2.7) 


( 2 . 8 ) 

(2.9) 
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We note G(x) is an even function, in the following, we only need to consider the 
behavior of G(x) when x —> 0 + . We divide our discussion into two parts. 

If s € (|, 1), we have 


1 f°° cos(x£) 

n Jo 1 + |£| 


1 


_ * f°° (sin ^) 2 \ 

2, d (=z(l 0 T+W- d( ~ 2 L TTlF d «) 


(sinf)^ 


2ssin(^) Wo 1 + |£| 
1 


2ssin(^) tt 
1 2 


9 r°° 

- -x 2 *- 1 / 




2ssin(^) tt 


^n- 


t 2s + x 2s 

t \2 


_ , _ 2s—1 | z „4s—1 

= do + aix H —x 


(sin|) 2 ar s (sin|) 2 \ 

~ps t 2s( t 2s +x 2sj) dt 

(sin |) 2 


tt J 0 t 2s (t 2s +x 2s ) 


d t. 


( 2 . 10 ) 


where 


Qo = 


TT \ ) Q 1 " 


(sin o 

-^—dt =-sr(—2s) sin(7rs) < 0. 

t 2s TT 


2ssin(§) 

We can write the last term on the right hand side of (12.101) as 


x 4s-1 (sin |) 2 / r x f 1 f°°\ a; 44 - 1 (sin |) 2 


fOO ^As — 1 Win t_\2 

_2 /_ Jy. _ 

/o t 28 (t 28 + a; 2s ) wo Jx 


a + r + f)i 


t 2s (t 2s +Z 28 ) 

For each term on the right hand side of (12.111) . we have 


dt. 


( 2 . 11 ) 


r ir 4s-1 (sin |) 2 
/o t 2s (t 2s + x 2s ) 


dt < 


Cx 2s 1 J ^ydt = 0(x 2 ), 


pl x 4s_1 (sin |) 2 


ni p 


4-ydt < Ca; 4s 1 / —j—dt = 0(x 2 + x 4s 4 ), 
t 2s (t 2s + x 2s ) - .4 t 4s v ; 


».4s—1 (t\2 

i o / o -W-di < Cx 4s , . 

t 2s (t 2s +x 2s ) Ji t 4 


X- "(sin |) 2 ^ <0 _ 4a _ 1 /' 00 (sin|) 2 


dt = O^ 48 " 1 ) 


Hence, we have 


x 48 Hsml) 2 ( min{2 , 45 _ 1}) 

3 t 2s (t 2s + X 2s ) 1 


as a result, we get the first one in (12.71) . 
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When s = |, we have 


cos(a’£) 




cos(^) 

l + £ 


d£ + 


f°° cos(:r£) 

* 1 + £ 


d£ 


_L_ 2(sinf ) 2 

■I + ? 1 + C 


d?- 


= 1 og( 1 + -)- 2 


= — log x + log — -2 


f 5 hhil! d( . 

o x+t 
r? (sin |) 2 


cos(a;£) 

TTF 

00 cost 


d£ 


d t 


t 


d t 


7 T X + t 
2 

pOO i 

/ cost 


t 


dt 


It is easy to see, 


+ 2x 


r 2 (sin |) 2 
o i(a; + i) 


dt —. 


(sin |) 2 


/ o + i) 


dt, 


cost 
t(x +1) 


cost 

t(a; +1) 


dt = 0(1). 


dt + O(x). 


Therefore, for s = we obtain 


G(x) = -logx + 02 + 0(x) as x —> 0 + , 

7r 


where 


02 = 


■ (log ■ 


- 2 


(sin |) 2 


dt 


cost 


dt 


( 2 . 12 ) 


□ 


3. The scheme of the proof 

Our strategy of the proof of the main results is based on the idea of solving the 
second equation in ed for v and then working with a nonlocal elliptic PDE rather 
than directly with the system. It is convenient to do this by replacing first u by 
r e u and v by t e v, which transforms ED into the problem 

( (—A) s u + u — ^ = 0 in R, 

< (— A) s v + e 2s v — r e u 2 = 0 in M, (3.1) 

D u,v>0, u,v—!> 0, as |a"|oo, 

with the choice of the parameter t e as in (11.101) . 



k 

u ~ u( x ~ 9<)» v ~ i> 
2—1 


Then, we have 
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i.e., the height of the bumps near the qi remains bounded as e —> 0. 

In the sequel, by T(h) we denote the unique solution of the equation 

(—A) s i> + e 2s v = r e h in R, v(x) -A 0 as |x| —y +oo, 

for h G L 2 (R), namely T = r £ ((—A) s + e 2s ) . Solving the second equation for 
in m we get v = T(u 2 ), which tends to the nonlocal equation 

„,2 


(—A ) s u + u — 


T{u 2 ) 


= 0. 


(3-3) 


We consider points q\, q 2 , ■ • • , qk, fe = 2 m in R which are the candidates for the 
location of spikes. These points (qi, < 72 > • • • ,qk) are restricted in the following set 


(91,92,-'• ,9fc) S A s = |( 91 , 92 , • • • ,9fc) I Qi = -Qk+i-i, 9i > 92 > • • ■ > 9fc, 


Let us write 


(,Ql1 Oli 1 Qm) ^ Qs,rj J* • 
k 

W(x) =Y,U( X ~qi). 


(3-4) 


We look for a solution to m in the form u = W + 0, where <f> is a lower order 
term. Then, formally, we have 

T(u 2 ) = T{W 2 ) + 2T(W(j>) + h.o.t., 

where h.o.t. corresponds to the higher order terms. We denote V = T(W 2 ). By 
G e (\y\) we denote the fundamental solution to (—A) s + e 2s in R. We can write 

k 

T(W 2 ) = r e f W 2 (y)G e (\x — y\)dy ~ r e ^ f U 2 {x - qi)G e (\x - y\)dy, 

i= 1 ■' 

where the integration extends over all R. For |x| = o(e _1 ) we have 

,1—2s I I 2s—1 I ~.3—2s\~\2 I _2s|,_|4s—1 


Ge{\x\) = 


otoE" ‘° + ai|x|" * + 0(e 3 2s \x\ 2 + £ 2s |x| 4s *), ifsG(i,l), 


-ilog( £ |x|)+0(l), 

Using this, and the definition of r e we get that near the qi s, 

V(x) = 1 + h.o.t.. 

Arguing similarly we get 

T(W(jj) = uj J W<j)+h.o.t., io = ^ . 

Then 

■ 2 W 2 + 2W(j>+ h.o.t. W 2 r 


if s = |. 


u 

v 


+ 2W())-2ujW 2 / W(j> +h.o.t.. 


V + 2T(W(/>) + h.o.t. V 
Substituting all this into (13.311 we obtain the equation for 4> 

(—A) s </> + (1 - 2W)(j> + 2 ujW 2 J W4> = S(W) + N{</>), (3.5) 

where S(W) = —(—A) S W — W + and N((j>) is defined by 

NW = 1WTW) ~ W - : 2W *- + 2 “"' 2 / w *’ 
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represents higher order terms in <j>. 

Thus we have reduced the problem of finding solutions to ( 13 . 11 ) to the problem 
of solving ( 13 . 51 ) for f>. We set = Zj. Rather than solving problem ( 13 . 51 ) . we 
consider first the following auxiliary problem: given points qj, find a function <f> 
such that for certain constants the following equation is satisfied 


Lcf> — S(W) + N((j>) + CjZj, (<j>, Zj) — 0, j — 1, 2, 


,k, 


(3.6) 


where 

L(j) = (—A) s 4> + (1 - 2 W)(j) + 2coW 2 J Wcf 

and (•, •) denotes the L 2 inner product. 

We will prove in section 4 that this problem is uniquely solvable within a class of 
small functions (j) for all points (qi, < 72 , ■ • • , qu) satisfying constraints (13.41) . Besides, 
the resulting constants Cj(gi, < 72 , ■ •' , ft) admit the expansion (11.111) . We will of 
course get a solution of the full problem whenever the points qt are adjusted in such 
a way that all of Cj’s vanish. We show the existence of such points in section 8, 
where the main result is finally established. In remainder of the paper we rigorously 
carry out the program outlined above. In particular, we will need to understand 
the invertibility properties of the linear operator L first. We will do it in the next 
section. 


4. The linear operator 

Proposition 4.1. Let U be the unique, positive, radially symmetric solution to 

{-A) S U + U -U 2 =0. (4.1) 

(a) There exists a positive constant d s depending on s only such that, as |ar| —>■ 
00, the following formula holds 


U(\x\) = 


|l+2s 


( 1 + 0 ( 1 )). 


Moreover, U'{x) < 0 for x > 0 and 


U\x) = - ^ + °( 1 )) as x -+ +00. 

(b) Let Lq = (—A) s + (1 — 2t/)id. Then we have 

Ker(L 0 ) = s P an {^}- 

(c) Let L be the operator defined in m and let 

L* = (-A ) s (j) + (1 - 2W)4> + 2 ujW J W 2 (j) 

be its formal disjoint. If we denote 

aw 

3 9qj ’ 
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then, for all j = 1 , 2, • • • , k, we have 

LZj = 0{ min | q, - qj^+^W (x) 2 + 0( min | Qi - qj \~^ s+ ^)W{x), 

*7*7 i j 

L*Zj = 0{ min \ qi - gi |-( 2 s + 2 ))W(:r). 

Proof. The proof of part (a) and (b) is given in [9] Proposition 1.1 and Theorem 2.3]. 
In fact, it was proven in [T] that (apart from translations) the function U(x) = 1+ 2 |i2 
is the only positive solution of m in H 2 (R) when s = 5 . Part (c) is a consequence 
of direct computations. □ 


We shall carry out the analysis of the linear operator L in a framework of 
weighted L°° spaces. For this purpose we consider the following norm for a function 
defined in R: given points q±, < 72 , • • ■ ,qk we define 

M* = ||p(aO"VlU~(R)> (4-2) 

where 

7 1 1 


p{x) = ■ 


< fJL < \ + 2s. 


We first consider a problem that will later give rise to the finite-dimensional 
reduction. Given a function h, ||/i||* < 00 find a f> and constants Cj , j = 1,2, • • ■ , k 
such that one has 

L(f = h + c jZj in R, 

4>{x) -A 0 as |x| —> 00 , (4-3) 

{(f,Zj) = 0 for j = 1,2, • • • ,k. 

By c we will denote a vector with components Cj. 

We refer to a pair (c t >, c) as a solution to (14.31) . We have the following existence 
result for (14.31) . 


Theorem 4.2. There exists positive numbers R and C such that, for any points 
qj satisfying \qi — qj\ > R for all i ^ j, and h with ||/i||* < 00 , problem i4-S\ ) has a 
unique solution = T(h) and c = c (h). Moreover, 

||T(/i)||* < G||/i||*. (4.4) 

The proof of Theorem 14.21 relies heavily on the following lemma. 

Lemma 4.3. Assume that qf,j = 1, 2, • • • , k, are such that min vti k" -Qj I -»• 00. 
\\hn ||* —> 0 and that <j> n solves 

L<j) n = h n + ^2 CjZj in 
j 

4> n (x) —t 0, as |x| —► 00 , 

(4>n, Zj) = 0 for j = 1,2, • ■ ■ , k. 

Then ||^ n ||* —>• 0. 

Proof. We will argue by contradiction. Without loss of generality we can assume 
that H^nll* = 1. Our first observation is that c" —> 0. Indeed, multiplying the 
equation by Zj and integrating by parts we get 

(<t>n, L*Zj) = c” f Z 2 + J2 c”< Zi, Zj) + (hn, Zj). 

J iAj 
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Using Proposition 14.11 by rather standard calculations, it follows that c" —> 0 as 
n —>• oo. 

Our next goal is to show 

J W<j> n —> 0 as n —>• oo. 

To this end consider the test function 

Z = -x ■ VW + 2W 
s 

and 

L 0 = (—A) s + (1 - 2U)id. 

We claim that 

LqZ = —2W + o(l). 

Indeed if we set u\(x) = X 2s W(Xx), then 

-(-A ) s u\(x) = \ 2s u\ -u 2 x + '^2 X is U(\x - &)U(\x - £j). 

*7 £ j 

Since Z = I|a=i the claim follows now from Proposition 14.II and the above. 
Decompose </>„ = a n W + ip n where ( W,ip n ) = 0. Then Lfi n = Loip n and we have 

o(l) = (L<t> n , Z) = a n (LW. , Z) + (LoV’n, Z). 


But 

and 


Mn, Z) = (ipn, L 0 Z) = -2(W, ifn) + o(l) = o(l) 


{LW,Z) = {W 2 ,Z) + o{ 1) = ^- J x,X7W 3 + 2 Jw 3 + o( 1) = (2-|-) Jw 3 + o( 1). 

It follows that a n —> 0, and ( W,fi n ) = o(l). Going back to the equation satisfied 
by <fi n , we see that 

(~ A ) s (f> n + (1 — 2W)4> n = o(l)(w 2 + ZjJ + h n = g n + h n , 

j 

with ||g ra ||* —> 0. By applying the Lemma 4.2 in [J, we get 
[1077.11* < C(\\g n \\* + ||h„||*) = o(l) 

which contradicts to the assumption ||<^>„||* = 1- Hence we finish the proof of this 
lemma. □ 


Remark: In the statement of Lemma 4.2 in U , it requires that g £ (|, l+2s). How¬ 
ever, we can go through the proof to get the same conclusion for the case g = l + 2s. 


Next we construct a solution to problem (14.31) . To do so, we consider the following 
auxiliary problem at first 

[ (~&) s <t> + 0 = g + £i=i CiZi, 

\ f Rl </>Zi = 0 for all i. { ’ 1 

Lemma 4.4. For each g with ||</||* < +oo, there exists a unique solution of problem 
<f> =: A[g\ S H 2s ( K). This solution satisfies 

ll-%][|* < C\\g\U. 
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For the proof of Lemma l4~4l one can see jU Lemma 4.3]. 


Proof of Theorem \4-2\ Let us solve problem (14.31) . Let y be the Banach space 

y ■■={<!>& ^(K 1 ) I II0b “ 11011* < °o}. (4-6) 

Let A be the operator defined in Lemma 14.41 Then we have a solution to (14.31) if 


we can solve 


We claim that 


<t> - A[2W<t>] + 2 A[ojW 2 j W<t>\ = A[g\, </> £ W (4.7) 


B[(j)\ := Bi[0] + B 2 [(t>] = A[ 2 W(t>] + A[- 2 ojW 2 J W<j>\ 

defines a compact operator in y. At first, we show that B 2 defines a compact 
operator in y. We begin to prove that C[<p\ = ujW 2 f W</> is a compact operator from 
y to y. Let us assume that 4> n is a bounded sequence in y. By direct computations, 
we get I(4> n ) = fW</>n is contained in a bounded set in M 1 . Therefore, we can pick 
out subsequence of /(</>„) converge, combined with the fact W{ x) < Cp(x), we get 
C is a compact operator from y to y. By Lemma |4.41 we have A is a continuous 
operator from y to y. Hence, we get B 2 is compact. 

Then, we come to consider the operator B\. We denote /„ = A[ 2 W(f> n ] and claim 
that 

I fn(x) - f n {y )| 


sup ■ 

Xy ty \ x V I 

By Green representation, we have 


< C\\<j) n 


= min{l, 2s}. 


(4.8) 


fn = A[2W<t> n ] = 2 J G(x, z)W<j> n dz, 


Since 


\\A[2W</> n ] || L c= = ||2 J G{x,z){W<f n ){z)dz\\ L ~ < C\\cj> n 
we can obtain 

I fn(x) - f n (y )| 


sup ■ 

x^y 


- y\ c 


< cWfa 


for \x — y\ > For \x — y\ < |, we have 

I fn{x) - fn(y )| < Cll^lloo f \G(x - z) - G(y - z)\dz. 

Jr 

Now, we decompose 

f \G(x ~ z) — G(y — z)\dz = f \G{z + y — x) — G{z)\dz 

Jr Jr 

| G(z + y — x) — G(z)\dz 


'|z|>3|y-x| 


'|z|<3|y-x| 


| G(z + y — x) — G(z)|dz. 


We have 

f \G(z + y — x) — G(z)\dz < f dt f \S7G(z + t(y — x))\dz\y — x\, 

J\z\>3\y—x | J 0 J \z\>3\v—x\ 


|z|>3|y— x\ 
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and since 3 |y — x\ < 1, 

f \VG(z + t{y-x))\dz<C(l + [ 

J\z\>3\v-x\ Jl 


dz 


12—2s 


< 


l>\z\>3\y—x\ FI 

<7(1 + \y - zl 28 " 1 ), « e (|,i), 

<7(1 + log \y — x|), s=|. 


On the other hand 

f \G{z + y — x) — G{z)\dz < 2 f 

J\z\<3\y—x\ J h 


< 


\z\<i\y-x\ 

C\x - y\, 


\G(z)\dz 


s e (311)> 


<7|x - y \log \x - y\, s = 


Hence, we get 


„ „ I fn(x) - fn(y) | / ^ MJ . 

T— 

Xyty y | 


Therefore, (Id.811 is proved. 

Using the Arzela-Ascoli theorem, we can get a subsequence of f n which we 
relabel the same, that converges uniformly on compact sets to a continuous function 
/. Outside such a compact set, by using the property of the Green function and 
standard potential analysis, we have 

f n = A[2W(j) n \ = ^ +2s (1 + o(l)) as |x| 00 , 

FI 

where 7 „ = 2y f R W<j > n . Based on the previous subsequence we have chosen, we can 
further choose a subsequence such that I{<j> n ) = / R W(f) n converges, combined with 
the fact p{x)~ 1 \x\~' 2s ~ 1 -A k as x —> 00 , we can get a subsequence of f n converges 
in y outside a compact set. This implies that B 1 is compact. Therefore we prove 
the claim that B defines a compact operator in y. 

Finally, a priori estimate tells us that for g = 0, equation (14.311 admits only the 
trivial solution. The desired result of Theorem 14.21 follows at once from Fredholm’s 
alternative. □ 


5. Basic estimates 

In this section we calculate basic estimates, including error estimates. Our first 
task is to analyze the solution V of the problem 

k 2 

(—A) S U + e 2s V = r e U(\x — qi\) , V(x) —> 0 as |x| -A 00 , 


2=1 

where r e is given in (13.21) . Denote by Zq the solution of 

(—A ) S Z 0 + e 2s Z 0 = U{\x\) 2 , Zq(x) -a 0 as \x\ -A 00 , 
and by % (x), » / j, that of 

(-A ) s 0 i:j +e 2s 0ij = U(\x - qi\)U(\x - q 0 1), 9ij(x) -A 0 as |x| -A 00 . 


(5.1) 
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Then we have 

k 

V(x) = T e ^ Z o(\% ~ Qi\) + Te ^ Oij(x). 

»=1 

We will now study Zq(|x|) in the range 

x G J s := [-¥£~,Me“] 
for s G ( 5 , 1 ) and in the range 

x G Ji := [—M(— loge)U M(— logs) 2 ] 
for s = i. Here M = —. 

Z T] 

By Green representation, we have Z 0 {\x\) = / G s {\x — y|)C/(|y|) 2 dj/. We can 
expand Zq{\x\) as 

7 ( \ ix / f (aoe 1 - 2 *+ a 1 \x-y\ 2s - 1 )U(\y\) 2 +0(e iain V+^- 2s ’i-^), s G (|, 

° \ -i/log(e|a;-y|)C/(|2/|) 2 + a2/C/ 2 + 0(e(-loge)^), s = 3. 

(5.2) 

Let us consider the quantity H s (x) — ai / \x — y | 2s 1 U(\y\) 2 for s G ( 5 ,1), by 
noting that U{\y\) is an even function, we can easily obtain 

H s {-x) = a l 1 1 - x - y\ 2s ~ l U{\y\) 2 = fll J \ - x + y^Udy]) 2 = H s (x). 
Furthermore, we can write 

H s (\x\) = ai|x| 2s_1 J U(\y \) 2 + / s (|x|), 

where f s and its first derivative are uniformly bounded. Similarly, we can get 

Hi(\x\) = -ik)g|x| J U{\y\) 2 + fi(\x\), 

where /1 and its first derivative are uniformly bounded. 

Let us now consider the function 0^ given in (EH). Since Oij can be represented 
as 

%=y G e {\x-y\)U{y-q i )U{y-q j ). 

At first, we consider the case for s G (^,1). For x G / s , the following expansion 
holds, 

Oij(x) =a ° £l 2S J u (v ~9i)U{y-qj) + a 1 J \x - y\ 2s ~ l U(y - q^U(y - qj) 

+ 0(e^) J U(y - qi)U(y - qj). 

Using Proposition 14.II one can show that 

J \ x ~ 2/| 2s— 1 £/(y — qi)U(y — qj) = 0(e^), 

uniformly in x G I s , a similar estimate holds for the derivative of the above expres¬ 
sion with respect to x. Let us set 

<S s (M) = J U(y)U(y - z)dy. (5.3) 
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Following a standard potential analysis, we have 

J,(N) = 0(l) N2fl+ 1 1 + 1 . (5.4) 

For the case s = \, we can expand the term Oij (x) as follows 

Qij(x) = ~~J lo s( £ l x - y\) u (y - <nW(y - <u) 

+ a 2 J U{y - qf)U{y - qj) + 0(e|loge|5). 

Using Proposition 14.II we can obtain that there is a S > 0 and close to 0 such that 

J {a 2 -\og\x-y\)U{y-q i )U{y-q j ) = 0( ^ J , 

uniformly for x £ 1 i, a similar estimate holds for the derivative of the above 
expression with respect to x. Let us set 

6i(\z\) = J U(y)U(y - z)dy. (5.5) 

Following a standard potential analysis, we have 

d<W) = (5-6) 

Thus, combing the above estimates we obtain: 

Lemma 5.1. For the term V, we have 

(a) When s £ (|, 1), the following estimate holds for x £ I s , 
k 

V{x) = l + T £ '^2H s (\x-q i \)+c s '^25 s (\q i -q j \) + 0(E 2s ~ 1 ). 

*= 1 

When s = the following estimate holds for x £ 1 1 , 
k 

V(x) = l + T e ^2Hi(\x- qi \) + ci X]^(|g^-gjl) + 0 ( _ 1 Q )■ 

i=l ijij 

Here c s ,ci are some generic constants independent of e 


A similar estimate holds for the derivatives ofV with respect to x. The 
function H s (\x\) is given in the above discussion and, for |x| > 1, has the 
expansion 


H(\x\) 


a iW 2s 1 fU{\y\) 2 + fs(\x\), s £ (0, |), 

-^°g\ x \f u (\y \) 2 + f±(\x\), s = \- 


The function S s is given in and \5.6\) for the cases s £ (i,l) and 

s = i respectively. 

(b) When s £ (^, 1) for a: 6 R \ I s , then we have the following lower estimate 


V{x) > C 


(5.7) 


1 + (e|a:|) 2s + 1 ' 







18 


JUNCHENG WEI AND WEN YANG 


When s = ~, for x £ ffi. \ Ii , then we have the following lower estimate 


V(x) > C- 


1 


1 4- (e|a;|) 2 ’ 

Estimate and (15.81) can be proven by using the potential analysis. 


(5.8) 


6. FURTHER ESTIMATES 


For brevity we shall denote Ui(x) = U(x — qf) and W = Our purpose 

in this section is to derive estimates for the quantity 

W 2 


S(W) = -{-A) s W-W + 

which can be written as 

W 2 k 

S(W) = — -Y, u l 

i= 1 

Our first result is the following 


V ’ 


( 6 . 1 ) 


Lemma 6.1. Let the number p. = 1 + 2s in the definition of || ■ ||*. 
qi satisfying conditions in and all sufficiently small e we have 

Ce s 4s , s £ (i, 1), 


For all points 


l|S(W)||*< 


s = i 2 ’ 

d 2 ’ 


where C is some constant independent of e and S is any small positive number 
independent of e. 


Proof. Let us assume first x € I s . We write 


S(W) = 


1-V 

V 


J2 U i+ 2V ~ X U i U 3 =Jl+J2- 


( 6 . 2 ) 


i =1 


^ 2 ’ 

s = I 
d 2 ’ 


At first, it is easy to see in the region under consideration, we have 

V(x) = | l + 

[ 1 + t A|log e | i -0> 

where S is any small positive number. Hence, 

' O(^XEli^), ^ e (1,1), 

0{ 1 


Jl = 


On the other hand, 

V-'UiUj < 2 UiUj < C 


■ I logs I 

1 


1 i^)(Eti^), * = £■ 

(i 


1 


I qi - 9j| 2s+1 M 1 + I* - qi |) 2s+1 
= OW--TT^i+T^®)' 


(l + \x-qj\) 


2s+l 


h - Qj 


J 2 


0{s s 4 1 , ) p {x ), se(|,i), 

0(p^)d(a;), s = i. 


Hence 
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As a conclusion, for x £ I s , we have 

\S(W)\ = 


0{e s *t)p(x), sG(i,l), 
°( |loge|i-* VW* s = i 


( 6 . 3 ) 


for small e. 

Outside the above region, let us consider the case s £ ( 5 ,1) first. Assume now 
a :€R\Jj. By using Lemma 15.11 we have 

k k 2s+l| |2s+l _i_ 1 

|s(ir)| < c( e a -« |x| 2 -+> +1) £ u ? < c £ A,) oW 

2=1 2=1 ' * W\/ 

= 0{e s ~^)p[x). 

For the case s = Assume x £ R \ /1 , we get 

k k 


(6.4) 


i SW |<c W + i)g(/?< cy: (i r ;w_ + J |)a 


p{x) 


= 0 {--^—)p(x) 
logs 

Combining relations ( 16 . 31 ) . ( 16 . 41 ) and ( 16 . 51 ) . we prove the lemma. 

Another quantity which will be crucial for the remaining arguments is 

1 = j S(W)Z i . 

We shall consider i = 1 only, since the other cases are similar. Observe that 
du^x-qx) _ _ du(x-qx) anc j ^hus we have 

oq 1 ox 


(6.5) 
□ 

( 6 . 6 ) 


,dU 


dU , 


-z = / a - w 1 E u ite( x - 91) + / ^ _1 E u * u ^ {x - 51) =h+h. 




dx 


We will estimate lx and I 2 separately. In fact, we will find the following expansions 


d 


h = ~ a s-Q^2U{q j -t?i)(l + o(l)), 


(6.7) 


& 1 


and 


h = 


f -0se a - 1 &Z# 1 to-9il a - 1 (l+o(l)), s£&l), 

\ -9 i l( 1 + °( 1 )), s = i 

Here a s , /3 S , s £ [|, 1) are some universal positive constants which are independent 
of e. 

We will establish (16.71) at first. Using Lemma T5. II we obtain 


dV x 


dUx 




&J 




dx 


Let us estimate f UiUj^- for i ^ j, we observe that if i,j ^ 1, then 

/ UiUj l£ = °^ lqi ~ qiU ~ 9l|)_ " (1+2S) )• 
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On the other hand, if i = 1, j ^ 1 we get 

For the right hand side of the above equality, by standard potential analysis, we 
can get that for a certain universal constant c s > 0 such that 

J U 2 {x - qi)U{x - qj) = t s U{q\ - qj){l + o(l)), 


with a similar estimate for its derivative. We leave the detail in the appendix. 
Hence, 


f V-^UiUj 



Qi) 


1 _d_ 

2 Cs d qi 


^2 u {q 1 -q j ){l+°(l)), 
i ¥-1 


and estimates (16.71) thus follows. 

Let us now consider the term I\. Following a similar procedure, we get 


h=~ 


/ k k 

U(x - ®) 2 { Y2 T eH s (\x - q 3 \) + CsY^ Ss(\qj - qi\)} 

i= 1 3 =1 




where H s and 5 S are given in Lemma 1-5. II Let us first estimate 

r d 

9iji= / U(x - qi) 2 S s {\qj - qi\) — U(x - gi)dx, 


with j l. For i = 1 this term is zero, while for i ^ 1 we can estimate, using (15.41) . 
(ESD and Lemma IO 

bud < C'bj - g i | _(2s+1) ( 1 + bi - gi|) _(2s+2) - 


Then we come to consider the terms 

lij = j U(x- q.i) 2 H s {\x - ~ 9i)- 

First we observe that the term corresponding to i = j = 1 vanishes, by symmetry. 
If i and j are different, and both different from 1, then the resulting term is of lower 
order, more precisely 

r g 

hj = U{x- qi) 2 H s (x - qj ) — U (x - qi) 

g r 

= - dqi J U ( x ~ 9i) H s (x - qj)U(x - q x ) 

f °( V >~5k(\ qi “%| 2s_lf/ fe -9i)), s€(i,l), 

\ 0( 1 )4r( 1 °gl% ~<lj\U(qi ~qi)), s = \. 

On the other hand, if i = 1 , 

hj = - J U{xfH s {\x - {qj - Qi)|)dm 

f b s -£-\q 1 -q j \ 2s - 1 (l + o(l)), sG(i,l), 

1 b |^r 1 °gb 1 -®l( 1 + °(i))> s = i 
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where b s are some generic constants and we used 

u( T \ = l a ° f U 2 (y)dy\x\ 2s ~ 1 (1 + o(l)), s<£(±,l), 

5 l -¥/ t 7 ’ 2 (y) d y lo sN( 1 + o ( 1 ))> s = b 

provided x is sufficiently large. We put the detail in the appendix. Now, as for In, 
we get 

hi = J U(x- (qi - qi)) 2 H s (\x\)^-(x) = 0(\q t - gi|“ 3 ). 

Combining the above estimates, we immediately get (16.81) . 

Hence we have found that 

Ej/i ^( a aU{qj -g 1 )+A»e 2a_1 |?j — h l 2s_1 )(1 + °(1))> s e (|, 1), 

[ a i U ^ ~Qi) + Pi l2S i!s M ) (! + °(1))> 

Thus, we obtain the following result: 

Lemma 6.2. For all points (qi,q 2 , ■ • • , qk) satisfies j3-4\) - If s £ (|, 1)) 


5(W)Zi = 


s = i. 




dqj 


-(1 + 0 ( 1 )), 


where 


ifs = b, 


where 


2s-lJ2s-l 


F s (r) = a s U{r) + /3 s e 2s V 


dFi (\qj - qi\) 

S(W)Zj = £ -^4^-(1 + 0(1)), 




dqj 


log r 

Fi(r) = aiU(r) + /3i- -p. 

2 2 2 log - 


7. The finite-dimensional reduction 

We will carry cut the finite-dimensional reduction process sketched in the first 
part of the paper. As in the previous section, we shall assume the points qt satisfy 
(IQ) . Recall from Section 3 that the original problem was cast in the form 

(-A) s u + u = ^ y (7.1) 

Rather than solving this directly we consider instead the problem of finding A such 
that for certain constants Ci one has 


(-A) S A + A 


(7.2) 
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and (A — W, Zi ) = 0 for all i. Rewriting A = W + <j> we get that this problem is 
equivalent to 


(-A) s (/> + (j) - 2Wf + 2 W 2 cj J W(j> 

W 2 (W + d>) 2 W 2 of 

= - (-A)' S W -W + — + - — - 2W<j> + 2 W 2 uj W0 + Y, dZi 

' i 

(7.3) 


V T{{W + 4>) 2 

= S{W) + N{<t>) + Y J CiZi 

and 

{<f >, Zi) = 0 for all i. (7-4) 

Using the operator T introduced in Theorem 14.21 we see that the problem is equiv¬ 
alent to finding a <f> G H so that 

<j> = T{S{W) + N{<j>)) = Q{<t>). 

We will show that this fixed point problem has a unique solution in a region of the 
form 


S, = 


I 


11011* < Ce s s G (|, 1), 


(7.5) 


v IWI.<C , T T5iip y ’ * = h J 

for any small positive constant 6, provided e is sufficiently small. Here 
H = {0GL°° | (<f>, Zj) = 0, j = 1,2, ■ ■ ■ ,k}. 

We recall that from Lemma EU 

/ \\S(W)\\. <Ce s ~^ : aG(|,l), 

\ IIW)IU<c n ^ 7 , s = l. 

On the other hand, N(<f)) admits the estimate provided by the following lemma. 
Lemma 7.1. Assume that cf> G V s . Then 

provided e is taken sufficiently small. Here 


a(e) = 


e s i, s <= (1,1), 
1 o — I 

1-61 ^ 0 7 


log e| 


as e —> 0. 


Proof. Let us assume first x G M \ I s and s G [2,1). We observe that using a 
standard potential analysis one can show that in this range of x we have W(x) < 
Cp(x) and 

T((W + « 2 ) > C 1+)i ‘ |)M , and T(W4,),T(4?) < C j - M|,. 

Then 

~2WVcj) + Vf 2 - 2W 2 T{W(j)) - W 2 T{cj) 2 ) 


M0)l < 


< C 


p{x) 2 


VT((W + (j)) 2 ) 
p(x) 2 


— 2Wcj) + 2uW / Wf 


+ 


.1 + (£| x |) 2s+1 1 + ( e | a :|) 2s+1 


■ 11011* 11011* + Cp(x) 2 \\<j)\\ 
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Therefore, we get 

\p(x)~ 1 N(cj))\ < C M*+o-(e) ll^ll*. (7.6) 

Let us consider now the case x £ I s for s £ [|, 1). We decompose N(<j>) in the 
form 

N((j>) = Ni(<j)) + N 2 {^), 

where 


Ni((f)) = {W + ^f 


1 


VT{{W + <j>) 2 ) V 


_ 1 + 2T(W<p) 


V 2 


- (2W + </>)</) 


2T(W(j)) 

V 2 


and 


N 2 (0) = -2</>W{l - i) + 2 W 2 [u J W(j>- 


T(W4>) , 4> 2 

v 2 J + V' 


We have that T((W + </>) 2 ) = V + 2T{W(j)) + T(</> 2 ). On the other hand, 

1 + 0(e s 4»), s 6 (f,l), 

1 + 0 (|log^|l-a)’ S = h 


V(x) = 


Also, 


T{W<j>) = 


to f W4>+ 0(e s_; fc)||(/)||*, se(i,l), 

u;/^ + 0( |log ^ |1 _i )||^||*, s=i 


and in particular |T(W<^)| = 0(||</>||*). Likewise, T(4> 2 ) = 0(||</>|| 2 ). Combining 
these facts we obtain 

|JVi(0)| < C(W + tfTtt 2 ) + C[2W$ + (t> 2 ]T(W(j>) 

<c P (x) U\\l 

A similar analysis yields, 

' Ce s -^{\<t>\W + W 2 1|0||*) + C|0| 2 , s G (§, 1), 

1 

2 ‘ 


\n 2 (4>)\ < 


I logs 


C t^MW + W 2 \\(t>\\*) + C\(j)\ 2 , s = t 


Hence, 




in this region. Combining this estimate with (17.61) . yields the result of the lemma. 

□ 


Using the definition of the corresponding norms, splitting different ranges of x 
as in the above proof, it is readily checked that the following holds: If 


II0* II* — 


Ce s ~~s , 

C(-uk) 1 - 5 ’ 


s £ (|, 1), 


i = 1,2, 


(7.7) 


then, given any small /i > 0, we can find e sufficiently small such that 

TOi)-iV(^)|U<H|01-^2||*. 


This implies that the operator Q is a contraction mapping in the set (D s defined 
in dm On the other hand, we also get from the Lemma 17.11 that Q maps T) s 
into itself. By using Banach fixed point theorem, we get the existence of a unique 
fixed point of Q in this domain, which depends continuously in the *— norm on the 
points of qi. We summarize this result in the following proposition: 
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Proposition 7.2. For all sufficiently small e and all points qi satisfying 1 3.4\ f 
we have the existence of a unique solution to 4> s = ■ ■ ,qk ) and c s = 

c s (qi , ,qk) which satisfies m- Besides, (4> s , c s ) depend continuously on the 
qi’s. 

In addition the following foi'mula holds for the components of c s j of c s : 


(7.8) 


- bs,j "1“ d 


s,ji j — 1) 2, • • • , 2m, 


with 



the error terms e SJ - satisfy 



Proof. We only need to prove the formula for c s j : s. Let us observe that the c s j 
satisfy the relations 



which define an ’’almost diagonal” system, from which the c s j’s can be solved 
uniquely. The main term in the above expansion is given by (S(W), Zf). To obtain 
the estimates for these numbers, which will equal the c s j’s at the leading order, we 
observe that 



and 



Formula (17.8|l is now an immediate corollary of Lemma 16.11 Lemma O and the 
expressions found for the c s j’ s. □ 


□ 


In the following section we will find the points qj such that all c 5;7 - ’s vanish, and 
satisfying the conditions in (13.41) . 


8. Solving the reduced problem 


In this section, we shall look for the point {qi,q 2 ,--- ,qk) such that c s j = 0 
and thereby prove Theorem 11.11 We first establish the presence of the zero of 
b s = (b s , i, & s , 2 ) • • • , b s ,k) and then use the degree theory to get the existence of the 
points (qi,q 2 , • • • , qk ) such that c s = 0. 

We recall that 



( 8 . 1 ) 
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where 

\ ot s U(r) + ^ s £ 2 s ~ 1 r 2s ~ l , se(i,l), 
r -W-{a,(/(r)+^, » = I 

It is not difficult to see that finding the zero point of b s is equivalent to finding the 
critical point of the following function, 

>9*0 = ( ot s U{qi - q 3 ) +/3 s e 2s_1 - 5 j| 2s_1 ) (8.2) 

for s £ (|, 1) and 

5i ( 51 ,( 72 , • • • ,5fc) = -9j) + /?ij—rl°gl% “ %l) (8-3) 

for s = i. 

Since qi and 5 fc +1 _j are symmetry with respect to the origin for i = 1,2, • • • , k. 
So, finding the critical point of S s , s £ [^,1) in (13.41) is reduced to finding the 
criticai point of S s introduced in (11.61) for s £ (^, 1) in (11.81) and (11.71) for s = in 
m respectively. 

For the functions 3 S , s £ [^, 1), we have the following property 

Lemma 8.1. The functions 3 S , s £ [^, 1) admit an interior minimal point in the 
set lw and m for s £ (^,1) and s = \ respectively provided p is sufficiently 
small. 

Proof. By Proposition l4.il we have as |a;| —> 00 , 

U ( x ) ^ | a ,|i+ 2 S ( 1 + °( 1 ))- 

For e sufficiently small, we have that S s admit the following asymptotic expan¬ 
sion 


S s (<7i,52,-- 

' , 5 m) — 

m U 

E d-d)) 

+ ^e 2fl - 1 |2g i I 25 ” 1 ) 




+ §(k-9fl 1+2s(1+ 

°(1)) + /3s£ 2s_1 |5i — 5i| 2s_1 ) 




y( b s a s 

^M9i+97l 1+2sl 

o(l)) + /3 s e 2s ~ 1 \q i + 5j| 2s_1 ) 

(8.4) 

for s £ (i, 1) 

and 




5i(9i,52, 

' ' ' i 5 m) 

™ /biQi 

(1+o(I))+ 

■^logi 10S|29il ) 




/ f>IQi 

(!)) + log 1 log 9j\) 




/ biai 

+ S(|« + «l> <1 + ° 

(!))+^i logl log|5i+5j|) 

(8.5) 


for s = 
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For convenience, we make the following substitution, 


\Qi~Qj \ = | 

Then, (18.4D - (18.5I) turns to be 

m 

Sa(qi,q2,--- ,9m) = b s a s £ S_ ^[^ ( 


£ 4 S S | di-dj\, SG(i,l), 


(log |) = |rfi -dj\, s = \. 


( 8 . 6 ) 


\2d ^ { l + °(l))+7.|2di| 2s ') 

+E L.-IiW 1 +°w) + * ~ d ^ 2s ~ 1 ) 


i= 1 


#3 


+ E(^-p^p - a ^ + °^)) + ls\ d i + d j\ 2s ')] (8.7) 


for s £ (i, 1) and 


1 m i 

,9m) = T [E(]2d-p(l + o(l))+7i log |2dj|) 

+ E (j5 ^f ( 1 + ° (1)) + logl * - dli ) 

+E L.+VpO-+°w)+^ lo s i*+^ i) 


+ (m 2 - -771)71 log(log-) 
2 2 £ 


(8.8) 


for s = Here j s = g-. 

Since the process of finding the interior global minimal point of S s for s 6 ( 1 , 1 ) 
in ciD and S s for s = ^ in m are the same, in the following we shall only give 
the detail of the case s = i. 

Before studying Si, we first consider the following function 
g(x’) = ic” 2 + 71 logic for x > 0. 

By analyzing the derivative of function g, we know that 

0( V = = k + ^ log2 E - pi lo §^' 

Let us come back to the function Hi. For (91, 92, • • • , q m ) € < 9 Q Sj7) , we have either 
there is some i such that di = 77 or di = i, or there are i,j such that |di — dj \ = rj. 
If the former case happens, i.e., there is some i such that di = rj or di = For 
convenience, we write 

Si(di,d 2 , • • • , d m ) = fbiai-E.'j 3i ((/i,92, • • • ,9m) - (m 2 - ^m) 7 i loglog 

2 V 2 2 jQgr ± / 2 2 2 •=■ 


: logH 
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Then, we can get 

5i (di,d 2 , ■ ■ ■ , dm) >min{i?r 2 (l+ o(l)) + 7 | log 2 r?, i?; 2 (l + o(l)) - 7 i log|} 

+ ( 2 m 2 - m - 1 ) [(^ + i log 2)71 - ^71 log 71 ] (l + o(l)) 
> Si (1,2, ■ • • , m) (8.9) 

provided 77 is small enough. If there are some i,j such that \di — dj\ = 77 . Then, 
Ei(d 1 ,d 2 ,--- ,d m ) > ( 2 m 2 - 777 -l)[(i + i log 2)71 - i 7 ilog 7 i] (l + o(l)) 

+ 77" 2 (1 + o(l)) +71 log 77 

> Si (1,2, • • • ,m) (8.10) 

provided 77 is small enough. By (18.91) and (18.101) . we get 

min Si (qi 1 q 2 , ■ • ■ ,q m ) > min Si (q 1 , q 2 , ■ ■ ■ , q m ). (8.11) 

qea<?i,, 2 qeg +n 

As a conclusion, Si (qi, q 2 , ■ ■ ■ ,q m ) admits a global interior minimal point in Q 1 . 

2 □ 

Before we give the proof of Theorem 11.11 we recall the following definition (see 
Definition 2.4 in [19] or in [?]). 

Definition 8.1. Let / : D —>• R be a C 1 — function, where D C R m is an open 
set. We say that Xq is stable critical point of / if V/(x 0 ) = 0 and there exists a 
neighborhood U of xq such that 

V/(x) ^ 0, Vz € dU, 

S 7 f(x) = 0 ,xeU +=+ f(x) = f{x 0), 

and 

deg(V/,t/,0)^0, 

where deg denotes the Brouwer degree. 

Remark: It is easy to see that, if xo is a global minimum point or a global maxi¬ 
mum point of the function /, then x$ is a stable critical point of /. 

Proof of Theorem \l.l\ By Lemma T 8 .11 we get the existence of (qi, q 2 , ■ ■ ■ ,qk) such 
that 

b s (qi,q 2 ,--- ,qk) = 0 . 

Furthermore, such a point is the global minimal point of the function S s in (11.81) 
and m for s £ (^,1) and s = \ respectively. We denote such point by q s = 
(?«,i, <3+2, • • • , q s ,m)- Now, we shall look for q = (<2T, q 2 , • • • , q m ) in the neighbor¬ 
hood of q s to make c s = 0. As we mentioned in the previous remark, such a global 
minimal point is a stable critical point of S s . Using the property of the stable 
critical point and relation between the function and 3 a , we can find an open 
neighborhood O qa of q s in Q StV such that the following holds 

deg(V qs 3 s ,O qa ,0) ^ 0 and V qa 3 s 0 on <90 qs , 
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which implies 

deg(b s , O qs , 0) ^ 0 and b s ^ 0 on <90 qs . 
According to the definition of , we can get 


|b s | > 


^ s+i- 

^ 2 2s 


s e ( 3 ; 1 )) 


> 1 - rr- . 

2 (logi)3 


on dO a 


2 ’ 


where s G [1, 1) are strictly positive constants. 
Next, Let us introduce the following homotopy, 

H(t, s, q) = b s + t{c s - b s ). 


We find that 

H(l, s, q) = c s and i?(0,s,q) = b s . 

It is known in Proposition 17.21 

( 0(l)e S+ ^~A, S 6(1,1), 7=;— 


Ic,o - bo I = 


m O q ,. 


o(l)(logi) 2 , s= 5 , 

As a result, we get H{t, s, q) ^ 0 on dOq s . Therefore, 

deg(iJ(l,s,q),O qs ,0) = deg(ff(0,s,q),O qs ,0) 


which implies 

deg (c s , O qs , 0) = deg (b. s , O qs , 0). 

We already know that the right hand side of the above equality is non-zero, there¬ 
fore, deg(c s , Oq s , 0) ^ 0. As a result, we can find q in O qs such that c s = 0. Hence, 
we finish the proof of Theorem ll.il □ 


The proof of even number bumps case is thus concluded. Let as assume that 
k = 2m + 1 and briefly sketch the way to proceed in this situation. In this case we 
introduce the set 

(21,02,-•• ,0/c) ^ K = {(01,02, -•• ,0/e) | 0/e = 0, 0* = -0fc_i, 01 > 02 > * * ' > 0/c-l, 

(01,02, , 0m) £ Qs^Y 

Now we consider the first approximation 

k 

W(x) =5>(*- ft ), (<7i, < 72 , ■ • • ,5fc) 6 A°. 

j=i 

In this case, because of the even symmetry, the ’’bad directions” corresponding to 
small eigenvalues are only those Zi with 1 < i < k — 1. With this observation made, 
the rest of the scheme of proof goes almost the same way. 
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9. Appendix 

In this section, we list the estimates used in previous sections and give a proof 
in the following lemma. 


Lemma 9.1. 


U 2 (\y\)U(\x - y\)dy = c s E/(|x|)(l + o(l)) as |x| oo, (9.1) 


J U 2 (\y\)\ x — y\ 2s 1 dy = c s |a;| 2s : (1 + o(l)) as |x|oo, 


(9.2) 


and 


J U 2 (\y\) log \x - y\dy = c s log |x|(l + o(l)) as |x| -> oo. (9.3) 

Here c s depends on the integral of f U 2 . 

Proof. Since the proof of (I9.1fl - (19.3f) are the same, we only give the proof of the first 
one. We divide the whole space into two parts, 

K = h U/ 2 := {y : \y\ < |x|3} u {y : \y\ > |x|t}. 

Then, for any y G Ii , we have 

1 11 1 / J/ y. 

I* - y\ 2s+1 = |cc| 2s+1 |1 - f P+ 1 = |a;| 2s+1 ^ + 

and 

=co + 0(\x\~^), 

where we used U(y) -> j^TTK- as \y\ oo and c 0 = f R1 U(y) 2 dy. Thus, 

[ U 2 (\y\)U(\x - y\)dy =U(x) [ U(y) 2 + f (U(x - y) - U{x))U{y) 2 
Jh Jh Jh 


h 

=c 0 L/(|cc|) + o(l)(|a;| _2s_1 ). 


(9.4) 


For y G I 2 , we have 


U{y) 2 U(x -y)< Cj^^U(x - y). 


Hence, 


f U(y) 2 U(x -y) = 0(1*1"^) = Od*!" 8 - 1 - 2 ^). 
■) h 


Combining (19.411 and (19.511 . we get EH)- 


(9.5) 

□ 
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